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Abstract 

We report on a lattice QCD calculation of the flavor-singlet axial coupling 
constant of the proton from the axial-vector current. The simulation is car¬ 
ried out at /3 = 6 on a quenched 16^ x 24 lattice. An extrapolation to the chiral 
limit shows that the connected insertion (valence and cloud parts) is 0.62±0.09, 
which is close to g\. The disconnected insertion (vacuum polarization from the 
sea quark) for the u or d quark is —0.12 ±0.01 and the s quark also contributes 
—0.12 ± 0.01. The total g\ is thus 0.25 ± 0.12 which is in agreement with 
experiments. In addition, we find g\ = 1.20 ± 0.11, g\ = 0.61 ± 0.13 and 
Fa/Da = 0.60 ± 0.02, also in agreement with experiments. 


PACS numbers: 12.38.Gc, 14.20.Dh, 13.60.Hb, 12.38.-t 



Recent experiments on polarized deep inelastic lepton-nucleon scattering from 
SMC [0 and E143 p| have confirmed the Ending of the earlier EMC results that the 
flavor-singlet g\ is small [Q. In the deep inelastic limit, the integral of the polarized 
structure function is related to the forward matrix elements of axial currents from the 
operator product expansion |^. Combined with the neutron and hyperon decays, the 
flavor-singlet axial coupling g\ is extracted. Since the axial current is the canonical 
spin operator, g\ is thus the quark spin content of the nucleon; i.e. g\ = Au+Ad+As, 
where the spin content Aq{q = u, d, s) is defined in the forward matrix element of the 
axial current, {ps\qi'j^'y 5 q\ps) = 2MNS^Aq. 

The fact that g\, which represents the quark spin contribution to the proton 
spin, is found to be much smaller than the expected value of unity from the non- 
relativistic quark model or 0.75 from the SU( 6 ) relation ( 3/5 of the isovector coupling 
g\ = 1.2574) has attracted a lot of attention. Despite numerous attempts to 
explain the smallness of g\ by means of various hadronic models, effective theories, 
and anomalous Ward identity 0, this problem has persistently defied satisfactory 
answer and has been dubbed the “proton spin crisis”. 

Attempts have also been made to calculate g\ via the anomalous Ward iden¬ 
tity using lattice QCD, which is well poised for such a challenge p, |^. However, it 
has been pointed out ||, that in the quenched approximation, the induced pseu¬ 
doscalar coupling can not be ignored due to the would-be Goldstone boson dominance 
in the forward matrix element. As a result, the topological coupling to the nucleon 
{N\GG\N) is not directly related to g\ in the quenched approximation. Lattice calcu¬ 
lations with dynamical fermions are free of this problem. Modulo uncertainties from 
the chiral and zero momentum extrapolation, a recent calculation of the topological 
coupling of the nucleon using dynamical staggered fermions 0 shows that the quark 
spin content is indeed small, comparable to the experimental findings. However, a 
calculation of g\ via the Ward identity does not address the issue why it is much 
smaller than that of the quark model, a question which is at the heart of the spin 
crisis. In an attempt to answer this question, we carry out a lattice calculation of the 
axial current directly. It turns out that the polarization due to the degrees of freedom 
outside of the naive quark model, i.e. cloud and sea quarks, are responsible for the 
smallness of g\. 


Lattice calculations of three-point functions have been used to study the EM 
axial (isovector) |^, and pseudoscalar( 7 riViV) [|^ form factors of the nucleon, 
the flavor-singlet g\, we calculate the following two- and three-point functions. 
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For 


G?.“p(i) = i:{0|x"Wr(0)|0) (1) 

X 

= E rf {o|x”(x/)A3(x)x'’(o)|o>, (2) 

Xf^X 

where is the proton interpolating field, Ls = —* 7375(1 -|- 74 )/ 2 , and A 3 is the 
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point-split axial current from the Wilson action 


Af, = 2nf{ma)[ip{x)^ij^'y5U^{x)ij{x + fi)+'ip{x + fi)^i'jf,j 5 Ul{x)'ip{x)]. (3) 


In contrast to the isovector case, the evaluation of the three-point function for 
the flavor-singlet current involves a disconnected insertion{D\) in addition to the 
connected insertion{Cl) 0]. The quark line skeleton diagrams for the Cl and DI are 
shown in Fig. 1(a) and 1(b) respectively. The DI refers only to the quark lines. They 
are nonetheless correlated via the background gauge helds. We shall calculate the Cl 
and DI separately. The factor f{ma) in eq. (^ is the hnite ma correction for the 
Wilson action. We take f{ma) = = {Akc/h — 3) for the connected insertion. For 

the disconnected insertion, the factor f{ma) is different from and depends on 
the Lorentz structure of the current. We shall use the value of f{ma) computed by 
comparing the triangle diagram in the lattice Wilson action and the continuum [|T^ . 


The connected insertion is calculated in the same way as the isovector axial cou¬ 
pling g\ where eqs. O and (^ are htted to two exponentials in the form 
iL jg-mq sixxiultaneously, using the data-covariance matrix to account 


fe-^^'f and 


for correlations. In evaluating the three-point function, the factor SKc{\TrUpiaqY^‘^ 
{kc = 0.1568) is divided from the axial current in eq. (|^) to account for the mean-held 
improvement of the lattice operator |]^ O, O. Numerical details are given in Ref. 


m- The unrenormalized lattice 5'^^con ^as been calculated for k = 0.154,0.152, and 
0.148, corresponding to quark masses of about 120, 200, and 370 MeV respectively 
(the scale a~^ = 1.74(10)GeV is set by the nucleon mass), and is plotted in Fig. 
2. The calculations were done on a quenched 16^ x 24 lattice at /3 = 6.0 with 24 
gauge conhgurations as in the isovector case O- Extrapolation to the chiral limit 
(kc = 0.1568) yields — 0.65 ± 0.09 as shown in Fig. 2. The g\ in the contin¬ 

uum is related to its lattice counterpart by the relation g\ = ZAgA", where Za is the 
hnite lattice renormalization constant. The one-loop calculation gives Za = 0.952 for 
/? = 6.0 IjT^, from which we hnd g\con = 0-62 ± 0.09. 

First we note that g\ = /S.Ucon+^dcon is the OZI preserving part of g\. If the DI 
part (sea-quark contribution) is roughly havor-independent, i.e. — As 

(strange quarks appear only in DI), as will be shown below, then gAcon should be 
close to g% the octet axial coupling, i.e. 


9a ^^con T ^d^Qn ~\~ (AUf^js T ^ddis 2As) — g a , 


( 4 ) 


From the recent ht of the nucleon and hyperon (3 decays, gA = 0.579 ± 0.021 
We see that our calculated gA con is quite in agreement with this. Second, we notice 
that the calculated ratio Ra = 9 a con/ 9 a is smaller than the expected value of 3/5 
from the SU(6) relation of the relativistic quark model. The behavior of this ratio 
Ra is plotted in Fig. 3 as a function of the quark mass. For heavy quarks (at ma = 
1, m = 1.74 GeV), Ra is 3/5. This is to be expected of the non-relativistic quarks 
where g^ = 5/3 and g/^ = ^- For light quarks, Ra becomes progressively less than 
3/5. It has been shown M] that this deviation from the SU(6) relation is due to 
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the presence of cloud quarks and antiquarks (those in the higher Fock space than the 
valence in the Cl). When this degree of freedom is eliminated by disallowing quarks 
from propagating backward in time, it is found that Ra becomes 3/5 (shown as dots 
in Fig. 3 where the errors are smaller than the dot size). The SU( 6 ) relation is 
therefore recovered in the valence approximation, in which the Fock space is limited 
to the valence quarks. We should mention that under this valence approximation, 
the quark spin content is still not equal to the non-relativistic value of unity. This 
has been explained [H] as a relativistic effect. Owing to the presence of the lower 
component in the Dirac spinor, a conhned quark with s-wave upper component still 
has a non-vanishing orbital angular momentum. Thus we conclude that the reason 
9A con is only about 60% of the non-relativistic value of 1 is because of relativistic 
effects and polarization due to cloud quarks and antiquarks. 

To calculate the DI in Fig. 1(b), we sum the ratio between eq. 
over t. It has been shown []^ that as tf » a, this sum becomes 


and eq. (i 


■zi: 

i=l t 




const + 


tf»a 


(5) 


Thus, we calculate the sum as a function of and take the slope to obtain the DI 
part of g\. Since the DI involves quark loops, and thus entails the calculation of off- 
diagonal traces of the inverse quark matrix, it poses a challenging numerical problem, 
for the size of the quark matrix is as large as 10® x 10® in the present case. In [Wl, we 
developed an efficient algorithm to estimate these traces stochastically with Z 4 noise, 
which gives the estimate with the minimum variance. This algorithm has been tested 
by computing quantities for which the exact answers are known. We have checked the 
cases of DI with vector and pseudoscalar currents. In both of these cases, the matrix 
elements are proportional to the 3-momentum transfer q. Hence, the forward matrix 
elements should vanish. Presented in Fig. 4 are results of the ratios in eq. (^ for 
the corresponding currents. These are obtained with 300 Z 4 noise vectors in each of 
the 50 gauge conhgurations for k = 0.148. To avoid contaminations by the boundary 
effects introduced by the hxed boundary condition we imposed in the time direction, 
we summed t from the nucleon source, which is 4 time slices away from the boundary, 
to 4 time slices from the other boundary. The slopes of the ratios (indicated as ME 
in Fig. 4) are htted from the point where the nucleon emerges as a single exponential 
in the two-point function which is at f = 8 and onward. We see that for both the 
vector and pseudoscalar cases, the calculated forward m.e. are indeed consistent with 
zero. We also calculated the scalar m.e. {N\qq\N)dis- We see that the slope is large 
and positive, and thus has the right magnitude and sign to remove the discrepancy 
between the vriVcr term calculated from ttN scattering and those lattice calculations 
with only the CL Similarly, the corresponding results for k = 0.154 are shown in Fig. 
4. This is the lightest quark we considered and it yields the largest errors. However, 
the forward m.e. for the vector and pseudoscalar currents are still consistent with 
zero and the scalar m.e. is larger than that of k = 0.148. We conclude from this study 
that stochastic estimation with Z 4 noise produces the correct results within errors. 
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Having tested our algorithm against known quantities, we proceed to calculate 
the axial-vector current. The results for k = 0.148,0.152 and 0.154 are presented in 
Fig. 5. They are obtained in the following way. First, the slopes are htted in the 
region tf > 8 where the nucleon is isolated from its excited states. The £t employs 
a data-covariant matrix to take into account the correlation among the time slices 
in the 50 gauge configurations and is htted over different ranges of to find the 
one with the minimum much in the same way the hadron masses are fitted. The 
per degree of freedom is small in each of the 3 cases considered here. They are 
given in Fig. 5. The resultant hts covering the ranges of tf with the minimum are 
plotted in Fig. 5. Finally, the errors on the £t, also shown in the hgure, are obtained 
by jackknihng the procedure. We see that the slopes are negative in the case of the 
axial current, in marked contrast to the scalar current in Fig. 4. 


In order to compare our results with experiments, we perform the extrapolation 
to the chiral limit. Plotted in Fig. 6 are the results of g\dis a single havor with 
the same sea-quark mass (denoted as Ki) as those of the valence- (and cloud-) quarks 
in the nucleon {K 2 ). These results include the one-loop renormalization constant 
as in the CL The DI has a two-loop log divergence in the continuum from the triangle 
diagram insertion on a quark line |^. We computed the hnite lattice renormalization 
associated with this 2 -loop contribution and found it to be much smaller than the 
one-loop result. Hence, we have neglected it here. The extrapolation is done with 
the covariant matrix to consider the correlation among the 3 k's. The error on the 
chiral limit result is again obtained by jackknihng the procedure of the extrapolation. 
To calculate As (the strangeness contribution to g\ we £x ki (sea-quark mass) at 
0.154 and extrapolate K 2 (valence-quark mass) to the chiral limit. These results are 
plotted in Fig. 7. 


From Fig. 6 and Fig. 7, we hnd that Audis = Addis = —0.12 ± 0.01 and = 
—0.12 ± 0.01. Together, we obtain g\dis = —0.36 ± 0.03. Combined with gAcow 'w® 
hnally obtain gA = 0.25±0.12 which is in good agreement with experiments We 
tabulate these and other results in Table 1 and compare with experiments. We hnd 
that they all are in good agreement with experiments. We should mention that our 
calculation on g^ is in agreement with a recent similar calculation with the volume 
source which predicts g]^ = 0.18(10). However, their g^, g% Fa, and Da are 
smaller than ours and the experiments by ~ 20%. This presumably is attributable 
to the fact that their lattice spacing at /3 = 5.7 is about 45% larger than ours at 
(3 = 6.0. We would like to point out an interesting observation. Comparing Fig. 6 
(ki = K 2 ) and Fig. 7 (ki fixed at 0.154, the strange quark mass), we notice that 
although the sea-quark mass (related to Ki) changes by a factor of 3 in Fig. 6 , the 
results still coincide with those in Fig. 7 for each of the valence-quark case. It shows 
that the DI depends sensitively on the valence quark mass but is independent of the 
sea-quark mass in the loop within errors. This is reminiscent of the hnding that, after 
the hnite ma correction for the Wilson action, the coupling between would-be U(l) 
Goldstone bosons is independent of the mass in the two 75 loops [Q for light quarks. 
The two currents are related to each other via the anomalous Ward identity and, 
in this context, the mass independence is consistent with a DI of 75 current being 
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dominated by the zero modes. It would certainly be interesting to try to verify this 
with an explicit calculation of the zero modes. 

Table 1: Axial coupling constants and quark spin contents of proton in comparison 
with experiments 


This Work Experiments 


g\ = Am -|- Ad -|- As 

0.25(12) 

0.22(10) i 

0.27(10)0 

g\ = Au — Ad 

1.20(10) 0 

1.2573(28) 


g\ = Au + Ad — 2 As 

0.61(13) 

0.579(25) 13 


Au 

0.79(11) 

0.80(6)i 

0.82(6)0 

Ad 

- 0.42(11) 

-0.46(6)i 

- 0.44(6) i 

As 

- 0.12(1) 

-0.12(4)0 

- 0.10(4) 0 

Fa = {Au — As)/2 

0.45(6) 

0.459(8) PI 


Da = {Au-2Ad + As)/2 

0.75(11) 

0.798(8) PI 


Fa/Da 

0.60(2) 

0.575(16) PI 



In summary, we have computed both the connected and disconnected insertions of 
the axial current in the proton in a quenched lattice calculation. Albeit they all agree 
with experiments, the systematic errors due to finite volume, discretization, renor¬ 
malization, and quenched approximation (which could be as large as 7% - 20%|T2 


will have to be addressed in the future. Nevertheless, the physical picture of g\ is get¬ 
ting clearer. The smallness of the quark spin content compared to the non-relativistic 
value of unity is, first of all, due to the fact that the combined relativistic effect and 
polarization of the cloud-quarks reduces the Cl to 0.62 ± 0.09, a value very close to 
g^, i.e. gjicon — 9a- This is because the DI is almost independent of the flavors u, d, 
and s. Furthermore, the sea-quark polarization is large and in the opposite direction 
of the proton spin. It is the sum of all these effects that produces a small g\. We 
should stress that our calculation is gauge-invariant in that no gauge hxing is applied. 
Since there is no gauge-invariant dimension-three axial operator for the non-abelian 
gauge held |^, , our result is not mixed with the gluon spin. Only the quark spin 

content contributes. In this paper, we have not attempted to determine the complete 
composition of the proton spin. That would entail calculations of the orbital angular 
momentum and the gluon spin, which we leave for a future investigation. 


This work is partially supported by DOE Grant DE-FG05-84ER40154. The au¬ 
thors wish to thank T. Draper, G. McNeile, A. Shapere, and G. Thron for helpful 
comments. 
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Figure Captions 

Fig. 1 (a) The connected insertion, (b) The disconnected insertion. 

Fig. 2 The lattice Qa^cou the connected insertion as a function of the quark mass 
ma. The chiral limit result is indicated by •. 

Fig. 3 The ratio Ra = 9a/9a tor the Cl is plotted as a function of ma. The results 
of the valence approximation are shown as •. 
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Fig. 4 The ratios in eq.(|^) are plotted for the vector, pseudoscalar, and scalar currents 
with quark masses k = 0.148 and 0.154. ME gives the htted slope. 

Fig. 5 The ratios of eq.(|^) for the axial current are plotted for the 3 k cases. ME is 
the htted slope. 

Fig. 6 The DI of the axial current as a function of ma. The quark masses in the 
valence and the sea are kept the same. The chiral limit result is indicated by •. 

Fig. 7 The same as in Fig. 6 except the sea quark mass (ki) is hxed at 0.154. 
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The connected insertion 



Fig. la 


for connected insertion 



ma=ln(4«:p//c —3) 
g-2 
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The disconnected insertion 



































